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Abstract 



The electric dipole moments of nucleons ( NEDM, ) are calculated 
using the method of QCD sum rules. Our calculations are based on the parity 
( P ) and time reversal ( T ) violating parameter in QCD and establish a 
functional dependence of the NEDM on 0, without assuming a perturbative 
expansion of this symmetry breaking parameter. The results obtained from 
the QCD sum rules approach are shown to be consistent with the general 
symmetry constraints on CP violations in QCD, including the necessity of: 
(1) finite quark masses, (2) spontaneous chiral symmetry breaking, and (3) the 
?7a(1) anomaly. Given the current experimental upper bound on the neutron 
electric dipole moment ( nEDM ), \d n \ < 10~ 25 e • cm, we find \9\ < 10~ 9 . 
This result is compatible with previous calculations of nEDM using different 
techniques and excludes the possibility of solving the strong CP problem 
within QCD via a dynamical suppression mechanism. 
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I. MOTIVATION 



In this paper, we study the electric dipole moments of nucleons ( NEDM, denoted as 
djsr ), which serves as an indicator of both parity ( P ) and time reversal ( T ) symmetry 
breaking |l[ . The main focus is on the possible violations of P and T ( or CP ) symmetries in 
the strong interactions, with Quantum Chromodynamics ( QCD ) as the underlying theory. 
In this picture, the nucleons ( N ) are treated as composite particles consisting of quarks ( q 
) and gluons ( G ) and the strong P and T violating interaction in the QCD Lagrangian is 
characterized by a 9 parameter ( see Sec.|l| for a definition of this parameter ). Such a 
P and T violating interaction among quarks and gluons generates a coupling of the nucleon 
spin to the external electric field, and the strength of this coupling is defined as the electric 
dipole moment of the nucleon. 

Our purpose is to establish a functional dependence of the NEDM on the 9 parameter, 
along with other fundamental parameters of the QCD Lagrangian, e.g., the quark masses 
m q , and the values of quark condensates R q . Based on the current experimental upper 
bound of the neutron electric dipole moment ( nEDM, denoted as d n ), which is less than 
10 _25 e • cm || , we can obtain an upper bound on the strong CP violating^ parameter 9. 
Previous calculations, based on effective models of QCD, require \9\ < 1CT 9 [Q. The puzzle 
of explaining such an unnaturally small number is referred to as the strong CP problem || . 

The problem is difficult and interesting, since an analytical calculation of low energy 
hadronic observables based on the QCD Lagrangian is not a trivial task. Furthermore, 
it turns out that an important property of QCD, namely - chiral symmetry, is closely 
related to the strong CP problem and imposes three stringent constraints on the possible 
breaking of P and T symmetries in QCD, which include: (1) the necessity of non-zero quark 
masses ||, (2) spontaneous chiral symmetry breaking f7|, and (3) the Ua{X) anomaly |§ 0. 
These symmetry constraints not only dictate the functional dependence of all CP violating 
observables on the QCD parameters, but also provide a dynamical suppression to the CP 
violating observables ||. 

A natural solution of the strong CP problem can be obtained without invoking a tiny 9 
parameter if the dynamical suppression is sufficient to diminish the CP violating observables 
below the experimental upper bound. Thus, we need to face the challenge how to realize 
these constraints explicitly in our calculations without employing a perturbative expansion 
on the 9 parameter. Our calculations on the NEDM problem, which is the first one based 
on the quark-gluon degrees of freedom and the QCD Lagrangian with a 9 parameter, will 
provide a critical answer to such an interesting scenario for the strong CP problem. 

This paper is organized as follows: An introduction to the NEDM problem and the 
strong CP violation is given in section O, where we also set up the notations used in this 
work. The hadronic and quark-gluon representations of the nucleon correlation functions ( 
NCF ) are discussed in section [TTI| and section [TV], respectively. The results obtained from 
both representations of the NCF are used to derive QCD sum rules, which are analyzed in 
section [V] and section |VT[ We conclude with a brief summary in section |V11| . 



Here we use the CPT theorem to translate time reversal noninvariance as CP violation. 
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II. INTRODUCTION 



Before we discuss the QCD sum rule calculations of the nucleon electric dipole moments 
1(|, it is useful to clarify some issues related to the strong CP violation in QCD. 

1. CP Violations in QCD 

Up to dimension four, the most general QCD Lagrangian in four dimensional space- 
time, consistent with Lorentz invariance, hermiticity, gauge invariance, is 

£qcd ^^iP^ + m$ e ie ^ + -G 2 + §^§GG (1) 

where P = + ig a B^ ) ■ Y (2) 

1 

and G^y = -e^ a pG aP , e i 2 3 = 1 (3) 

The meanings of various symbols are: 
ip : quark field 

ijj : Dirac adjoint of the quark field, ip = V^7o 
B^ : gluon field, a = 1, .., 8 
A a 

— : generators of the color SU(3) gauge group, a — 1, .., 8 
Gf, u : gluonic tensor field, G^ = [ <9 M + ig 8 B^ d v + ig 8 B v ], G 2 = G^G^ 



9s 

0G 



strong coupling constant in QCD 
quark chiral phase 
gluon chiral phase 



Here we have two P and CP violating ( but C even ) terms in the QCD Lagrangian, 

im q sin 9 q ^TsVs -^T1 GG 

the former is referred to as a quark pseudo-mass term, the second is referred to as a 
gluon anomaly term. Our notation is chosen such that the corresponding CP violating 
parameters 9 q and 6q are angular variables^ and the QCD generating functional is 
periodic with respect to these CP violating parameters. The general QCD Lagrangians 
with two chiral phases 6 q , 6q are not all physically independent. Through the {7a (1) 
anomaly [11] [12|, we can shift some part of 9 q to 9q, and vice versa, by performing 



an {7a (1) rotation on the quark field 



2 The gluon anomaly term, when evaluated with instanton configurations, gives integer values. 
Thus, the QCD generating functional is a periodic function of 9q. 
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if) -> V' = e^ 75 ^, ^ ,t = -0t e «T» j (4) 

^ - g - 20, G - G - 20. (5) 

Therefore, only the difference between these two phases, 

= e G - e q (6) 



which is invariant under an Ua(1) rotation, is a physical parameter and can be used 
to label the equivalent classes of CP violating QCD Lagrangians ||. Furthermore, 
since physical observables should be independent of the reparameterization of the 
Lagrangian, we conclude that a CP violating observable should only be proportional 
to 9, instead of being an arbitrary function of 9 q and/or 9q. 

In addition to the explicit symmetry-breaking parameter 9, there are three important 
symmetry constraints which could suppress the magnitude of a CP violating observ- 
able. These constraints include: (1) explicit chiral symmetry breaking due to finite 
current quark masses, (2) spontaneous chiral symmetry breaking with a nonzero quark 
condensate, and (3) the Ua(X) anomaly ||. Suppression of the CP violating observables 
through the symmetry constraints is possible because in any of these particular limits, 
the effects of strong CP violation vanish, even with a nonzero 9 parameter. Therefore, 
it is desirable to have a calculation of a CP violating observable, e.g., d N , without as- 
suming a perturbative expansion in the 9 parameter. Such a calculation could provide 
an answer to the important question: Why is strong CP violation small? 

2. Nucleon Electric Dipole Moment as an EM form factor 

In a nondegenerate system like the neutron, the existence of an EDM implies the 
violation of both parity ( P ) and time reversal ( T, or CP ) symmetries. To establish 
a connection between the CP violating parameter 9 in QCD and the NEDM, it is 
useful to study the nucleon EM matrix element: 

V?(q; Pl ,p 2 ) = J d^x e**» (NbdlJuWNfa)) 

= (2vr) 2 ^(g-p) (N(p 2 )\.U0)\N( Pl )) 
= {2iz) 2 «5 2 «{q- V )V»{p 2 , Pl ) (7) 
with p = p 2 — pi = q (8) 

and 2uj is the space-time dimension. The nucleon EM vertex is extracted from the 
nucleon EM matrix element by factoring out the ( on-shell ) Dirac spinor for the 
nucleon state: 

v?{ P 2,Pi) = u( P2 ) r> 2 ,pi) «(pi) (9) 

Using (1) current conservation d^J^ 1 = and (2) hermiticity V^(p 2 ,pi) = 
[ (P2,Pi) ] , we can write down a general form for the EM vertex r^(p2,Pi) of 
spin 1/2 on-shell nucleon state: 
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n 3 f / ! '-M/v 

(10) 



r> 2 ,Px) = W) 7 , + iFf (g 2 )|^ - F 3 "(g 2 )^p + (g 2 )(g 2 7 , - g M g) 7s 



where we have four form factors , F^, F^ , and i 7 ^ characterizing the EM properties 
of the nucleons. At q 2 = 0, they are the various EM moments^ of the nucleon state: 

e F»{q 2 = 0) = Q N ( charge ) (11) 
[F^g 2 = 0) + Ff (g 2 = 0)] = fi N ( magnetic moment ) (12) 



2M N 

[F^(g 2 = 0)] = (In ( electric dipole moment ) (13) 



2M, 



N 



[F^(q 2 = 0)] = a^v ( anapole moment ) (14) 



2M 2 N 



It is useful to notice that the tensor structure associated with the anomalous magnetic 
moment F^ ( ) and that associated with the electric dipole moment F^ only 
differs by a factor z 75 . In view of this, we can rewrite these EM form factors in a polar 
form, 

F 2 + = F N e iaN ^ (15) 

(F N ) 2 ee (F 2 N f + (F 3 N f (16) 

taiaa N = -^TT (17) 



Calculation of Hadronic Matrix Elements from QCD 

At first sight, to calculate nucleon EM moments from QCD requires a knowledge of 
the nucleon wave function in terms of the quark gluon basis and a technique of solving 
the nontrivial strong coupling dynamics in the low energy region. This is certainly 
beyond our current ability ( except for numerical lattice calculations ) and we have to 
rely on other approaches to avoid these complications. For this purpose, we choose to 
calculate a nucleon correlation function ( NCF ) in the presence of external EM fields. 

U N (p) = J d 2 "x e*»< T VN (x) fj N (0) ) eqfiG , F , v (18) 

where the nucleon interpolating field r/jv is a composite quark operator carrying the 
same quantum number as a nucleon. For our calculation, we choose 

n n = (d*cv) 75 7^ (19) 



3 Due to the use of a constant EM background, the contribution of the anapole moment F± 
vanishes in our calculations. 
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as a neutron interpolating field |13] . A similar expression for the proton can be obtained 
by exchanging u and d quarks. 

In this approach, the EM matrix element of nucleons, which describes the response 
of the nucleon states to the weak external perturbation, can be imbedded in the first 
order expansion of the NCF with respect to the external field. Since we shall focus 
on the EM form factors at zero momentum transfer, the EM fields can be taken as 



constants 14 



To the first order of the electric charge, the NCF Hn{p) can be expanded as 

U N (p) = 11$ (p) + eU^(p)F, u + 0(e 2 ) (20) 

We shall call \p) a nucleon propagator and IT^(p) a ( nucleon ) polarization tensor. 
The former describes the propagation of hadronic states carrying the nucleon quantum 
numbers, the latter gives the EM vertices of hadronic states with the external fields. 

4. Tensor Structures of the Nucleon Correlation Function 

It is important to know how to write down a complete set of covariant tensors ( these 
will be referred to as basis tensors, composed of one Lorentz vector p M and 16 Dirac 
matrices ) and decompose the NCF in terms of these basis tensors. Such structures 



come out naturally from both hadronic representations ( see the discussions in Sec. [Ill 
) and QCD calculations ( see the discussions in Sec. |TV| ) for the NCF. Furthermore, 
the QCD sum rules are extracted from the coefficient functions associated with these 
basis tensors. A correct decomposition will help assure that there are no omissions 
and redundancies in our calculations. 

We find that it is convenient to use a commutation-anticommutation relation analysis 



[fyl to generate all possible independent invariant tensors for the NCF. For the nucleon 
propagator 11$ (p), there are 4 independent tensors: 

I, 7s, P = P^- 7 M , P75- (21) 

For the polarization tensor Yl^^p), the tensor basis consists of 8 independent second 
rank tensor matrices: 

a^, a^-7 5 ( 22) 

PY-PY, (pV-pV)"75 (23) 

e^ Pa lp, e^p alp ■ 75 (24) 

- pV), Pip^r - p u r) ■ 75 (25) 

Given these classifications of the basis tensors for the NCFs, we can define various 
coefficient functions associated with them. 

U N (p) ee nS?(p) + eIl£»iV + 0(e 2 ) (26) 

ni?(p) = JT(P 2 ) ■ P + f 2 N ( P 2 ) ■ I + ^f 3 N (p 2 ) ■ 75 (27) 

= fi(p 2 )-P + fitf)-e»«<&» (28) 
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with 



lf2 N (p 2 )? 

tan^Ar(p 2 ) 



if?(p 2 )} + tir {p 2 )\ 

h N ( P 2 ) 



f 2 N (p 2 ) 



(29) 
(30) 



n£(p) 



+ #f (P 2 ) v{P^l v -P v l^) 



//j v (/rl <r'" r 



+ #f(p 2 ) p(p M 7^ - p v ^)e itp " (p2)75 



+ ^V) 

+ 3f(p 2 ) ^(^7^-^)75 
+ <7eV) #(p"7" - P^hs 



(31) 



with 



[sfV)] 2 
tan^Qo 2 ) 

[^(p 2 )f 
tan^(j9 2 ) 



[fitfV)] 2 + [^(P 2 )] 2 
[^(p 2 )f + t^(P 2 )] 2 



(32) 
(33) 
(34) 
(35) 



Notice that because of charge conjugation symmetry, the coefficient function associated 
with the P75 tensor in the nucleon propagator and those associated with the p^^f — p v ^ 
and e^^PcTfp tensors are identically zero. 

The discussion in the next two sections will be devoted to the constructions of a 
hadronic and a quark-gluon parameterizations for all these invariant coefficients. 

5. CP Violating Vacuum Condensates and the Quark-Gluon Chiral Phases 

In the QCD sum rule []TB[ calculations of the hadron observables, we use an operator 
product expansion ( OPE ) to expand a hadronic correlation function in the order of 
operator dimensions. The nonperturbative ( to be precise, nonanalytic in the strong 
coupling constant ) structure of the QCD dynamics is parameterized in terms of various 
vacuum condensates and the perturbative contributions ( Wilson coefficients ) can be 
calculated using Feynman rules. By truncating the OPE series at a certain dimension, 
we obtain an approximate representation of the hadronic correlation function in terms 
of QCD parameters. On the other hand, the hadronic observables can be built in the 
hadronic correlation function by inserting a complete hadronic states and expanding 
the hadronic correlation function according to the hadron invariant masses. Through a 
matching between quark-gluon and hadron representations, the values of ground state 
observables can be extracted. It is crucial that we obtain the values of various vacuum 
condensates from other sources rather than calculating them within the sum rule 
method. This becomes a problem if we need to include higher dimensional condensates, 
which are poorly known. Also, in the external field method, due to the polarization 
of the QCD vacuum, there appear so-called induced condensates we need to take into 
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account [[H|. If CP is not a good symmetry of QCD, there could be in principle more 
unknown condensates associated with the CP violating operators, e.g., i{q^^q) dq Bq and 

(GG) e q g . However, it is possible to relate these parity doublet condensates ( {qq) e Bq 
and i{q^^q) e 9a ) to the chiral phases we have introduced. Indeed, there is a simple 
theorem, which can be used to relate quark condensates in the parity doublets: 



Theorem 1 



(qq)t 



ifia 



+ 



i(qi5q)t 



„8 G 



is invariant under Z7a(1) chiral rotations. 



The above theorem implies that the two real number (qq) e Bc and i(q r Ysq)$ g G can be 
thought of as the coordinates of a two dimensional plane, with R q defining the radius 
of a chiral circle generated by the Ua(\) chiral rotations. 

(qq)e q ,e G = -R q cos 9 G (36) 
i(qi5q) dqt e G = -R q sm9 G (37) 

For this reason, we shall refer R q as a chiral radiusQQ. Further use of this theorem will 
be discussed in latter section, see Sec.[TV|. 

6. Ua(1) Chiral Rotations and the Use of the Polar Form in the Sum Rule 
Calculation 

We have seen in the previous discussions that we can rewrite many variables and/or 
parameters in the polar form. These include: 

(a) the quark ( 6 q ) and gluon ( Qq ) chiral phases in the general CP violating QCD 
Lagrangian ( Eq.(P ), 

(b) the anomalous magnetic moment ( ) and the electric dipole moment ( ) 
of the nucleon ( Eq.([15|) ), 



4 The value of a chiral radius depends on the quark mass. For light flavor m q < Aqcd, spontaneous 
chiral symmetry breaking implies R q is a finite positive number in the massless limit. 

5 For more than one light quark flavor, the values of the chiral phase 

Oqo = arctan ^ B ifiG 
w)o q fi G 

can be determined from Crewther's condition IliT 



$>G« = 0G (38) 



q 



m qi R qi sm(6 qi - 9 Gqi ) = m qj R qj sm(6 qj - 9 Gqj ) (39) 
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(c) the invariant coefficient functions associated with the tensor basis ( 

(d) the quark condensates in the parity doublet ( Eqs.(|36|),(|37|) ). 

We find it is quite convenient to adopt this convention for the following reasons: 

(a) The reparameterization invariance of physical observables under the Ua(1) chi- 
ral transformations of the QCD Lagrangian can be maintained throughout our 
calculation. Thus, there is no need to stick to a particular representation of the 
QCD Lagrangian. 

(b) The polar form allows natural identifications among chirally invariant and/or 
chirally covariant variables. Those chirally invariant observables depend only 
on the chirally invariant parameters, e.g., m q , R q , and 9; the chirally covariant 
variables change by a constant phase under Ujy(l) chiral rotations. 

(c) Since the 9 parameter appears as an angular variable, one can solve the 9 depen- 
dence of the CP violating observables without using a perturbative expansion. 
Furthermore, the periodic structure of the 9 dependence comes out automatically 
due to the polar form. 

(d) The symmetry constraints on the strong CP violations in QCD can be made 
transparent in the sum rule relations if we organize both the OPE series and the 
hadronic representations in the polar forms. Without an explicit solution to the 
sum rule relations, one can show that CP violating observables vanish if chiral 
symmetry is exact in QCD. 



III. STUDY OF THE NUCLEON CORRELATION FUNCTION ( NCF ) FROM 
THE HADRON DEGREES OF FREEDOM 



1. Basics 

One advantage of the QCD sum rule calculation is to extract ground state observables 
from a correlation function without knowing the exact wave function of the nucleon 
state. The price for this convenience is that the interpolating field we choose couples 
to all possible hadronic excited states with the same quantum number as the nucleon. 
In addition, the ground state matrix element we are interested in is often accompanied 
with other excited state contributions to the correlation function. Consequently, the 
extraction of ground state observables from a NCF is possible only if we can identify 
and isolate various contributions to the NCF from hadronic states. To achieve this 
purpose, we insert a complete set of hadronic states in the NCF Hn{p) of the inter- 
polating field 7/tv- In doing so, we can factorize the correlation function into nucleon 
spinors and hadronic matrix elements. While the ground state observables, e.g., mass 
M/v and EM moments F t N of a nucleon N, can be specified explicitly, the fine details of 
the excited state spectrum are smeared out by employing suitable parameterizations 



14]. The ground state nucleon observables, together with the excited state parameters 
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are basic ingredients of a hadronic representation of the NCF. We shall discuss the con- 
structions of hadronic representations of the nucleon propagator and the polarization 
tensor in the following two subsections. 

Since we wish to maintain an {7a(1) reparametrization covariance in our calculations, 
it is important to keep all chiral phases explicit. In particular, we fix the chiral phases 
of the physical hadron states ( and the QCD vacuum, denoted as ft ) to be zero and 
allow quark fields to be in any chiral basis. With this point in mind, the definition of 
the nucleon spinor is given by 



rj N | N(p, s N ) ) dq 0G = X N e l — 75 u(p, s N ) (40) 



Here the nucleon residue A at gives the overlap between the nucleon state N and the 
interpolating field rj^, the nucleon chiral phase 9n specifies the quark basis of the QCD 
Lagrangian. It is convenient to choose an interpolating field which transforms like a 
quark field under an Ua,(1) rotation. For instance, for our chosen neutron interpolating 
field rj n ( Eq.pD ) 



u — > u - l ^ 75 



d^d' = e^d (41) 
Vn^v' n = e l ^ 5 Vn, 6 n - 9' n = 9 n + 9 U (42) 

In this case, the nucleon chiral phase On transforms covariantly ( changing by a con- 
stant phase ) and the nucleon residue \n stays invariant under an £7a(1) rotation. 

2. Hadronic Representation of the Nucleon Propagator 

For the nucleon propagator EE^ (p), we can insert a complete set of hadronic states 
J2n 1^0(^1 = 1 between the time-ordered product of the interpolating fields r}N,f}N- 

UP( P ) = J (n\T( m (x),fj N (0))\n) eq>eG (43) 

= £ / d^x e *» e(x )(n\ VN (x)\N) eq!e jN\fj N (o)\n) eq>9G 

N J 

- e{-x Q ){n\n N mN) eqfiG {N\ VN {x)\n) eqfiG (44) 

With the definition of the nucleon spinor in a general quark chiral basis ( Eq . (|40|) ) 
and the standard procedure to simplify the algebra, we obtain a propagator of the 
nucleon state, with an overall chiral conjugation^ 

6 P-M N 6 ~ Xn { f-M* ) (45) 



6 A chiral conjugation of a two-point Green's function is defined as G(p) — > U{9) ■ G{p) ■ U(6), 
with U{6) = e^ 75 . 
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Similar expressions can be written for excited states N* with different residues \n* 
and total masses M^* . 

The contributions of the excited states to the NCF lead to many unknowns in our 
calculations^. Consequently, we shall take a simple parameterization to replace the 
contributions of all excited states without involving the complete hadronic spectrum. 
Based on a duality argument, we can identify the total contributions of the excited 
states to the nucleon propagator as the the leading terms from the quark- gluon calcu- 
lations, starting from a continuum threshold Sq : 

^ „/v/ 9n / \ 1 „ f 00 Ira f, N (s) ( quark-gluon ) , . , 

Re f? (p ) ( continuum ) = - Pr. I Jl V ; y ds (46) 

7T J Sq p S 

Here we use a dispersion relation for the invariant coefficient functions ( See 
Eqs.(p7|),(E28|),(|3lD ) to relate the contributions from two representations and Pr. means 
principal value of the complex integral. The calculations of the quark-gluon represen- 
tation of the NCF will be given in the next section. 

The hadronic representation of the nucleon propagator U$ (p) is given by the sum of 
Eq.(f45|) and Eq. (|46|) , where we have three nucleon variables, \^,9 N , and M N to be 
determined from the QCD sum rules. 

3. Hadronic Representation of the Polarization Tensor 

Since the polarization tensor n^T(p) comes from a time-ordered product of the electro- 
magnetic current J M and the interpolating fields rj N ,f] N , we insert two complete sets 
of hadron states, £tv \N)(N\ = £iV' \N')(N'\ = 1. 

eU^(p)F^ = J d 2 "x d 2 "z ( T£ int (z)VN(x)fj N (0) ) dqfia (47) 

= E E / d2 " x eWX I d2 " z d ( x - z ) d ( z - °) x 

N N' 

+ time ordering (48) 
where the vector potential N 1 of a constant EM field is given by 

A a ee -xpF Pa (49) 

The two insertions are independent. Hence, the polarization tensor H^{p) contains, in 
addition to the (1) nucleon ground state EM form factors, (2) ground state to excited 



This problem can be traced back to the choice of an interpolating field for the nucleon state. A 
perfect choice of an interpolating field is equivalent to the solution of an exact wave function of the 
nucleon state, which couples only to the nucleon state and has zero overlap with any excited state. 
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states transitions, and (3) transitions among excited statesQ. We can think of the 
polarization tensor, in its hadronic representation, as a huge matrix in the Fock space 
and the three contributions are given in a tabular form ( see Table 1 ): 

We shall discuss these contributions to various invariant coefficient functions of the 
polarization tensor, with the superscript labels g^h 9m^ 9m ( * * s the index for the 
independent tensor basis ) corresponding to three regions listed above. Thus, 

9? = 9m + 9m + 9%l 2 = 1, 2, 6 (50) 



the nucleon ground state EM form factors 

If we only take the ground state nucleon from the double sum over hadronic 
complete sets, the polarization tensor reduces to a product of two ( CP conserving 
) nucleon propagators and the nucleon EM vertex, with an overall chiral 
conjugation. 



e n^ 1/ (p)( nucleon state ) F M „ 



dTz e l * 
+ time ordering 



P- Af 



A? 



p-M N 



e* 2 



(51) 



Substituting all ingredients ( Eqs. fllOD , fllCf) , ([49]) ) into the equation Eq . (|5T|) , we 
obtain the contribution from the ground state nucleon to the polarization tensor 
( We have factored out an overall constant 



4M N (p 2 



g$\(p) = 2icos6 N M 2 N F l N + i CO s6 N (p 2 + M 2 N )F? - sm6 N (p 2 - M 2 N ) 
g^lip) = 2% sin 9 N M 2 N F? + i sin N {p 2 + M 2 N )F? + cos d N {p 2 - M 2 N ) 



f) 



g W(p) = 2i sin N M 2 N F* + i sin 6 N (p 2 + M 2 N )F? + cos 9 N (p 2 - M 
<#>(p) = 2M N (Ff 
g^lip) = 2M N F 

g%l(p) = (-2) [ cos 9 N F^ + sin 9 N F, 



)F» 



T + F, 



9m 



(p) = (-2) [ sinM^ + cos^K 



?N 
3 

.AT 
3 



= (-2) 
= (-2) 



F N cos( 
F/v sin( 



i N - a N ) 
V - oi N ) 



(52) 
(53) 
(54) 
(55) 
(56) 
(57) 



(b) ground state transitions to excited states 

As in the case of excited state contributions to the nucleon propagator, we need 
to sum over all ground state to excited states transitions, and parameterize the 
spectral function with a few constants. This can be achieved if we take the 
following expression: 



8 The transitions in a constant EM background are possible because we are looking at highly 
off-shell states. 
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( a 1 7^(0) | n ) 



n | 77iv(0) | N ) f 



= / d 2 »z ( O | T (£ nt (z) ^jv(O)) | N(p,s N ) ) t 



fia 



,0 G ,F^=0 

— pole term 



= eF 



+ 



— g J 2 < 5 



p^ u - p u 7» 



A' 



(58) 
(59) 



This is the most general form for a nucleoli spinor in the presence of an external 
EM field; we use the Dirac equation for the ( on shell ) nucleon spinor u(p, s^) 
to reduce the invariant tensor structure with four unknown model parameters: 
E^,(p N ,E^,(p N . All of these parameters are invariant functions of Q 2 = —p 2 . 
The independent tensors are chosen such that Eft, E^ are invariant under an 
Ua{1) rotation, and the phases ip N , <p N transform covariantly. As a crude approx- 
imation, we shall neglect the Q 2 dependence of these model parameters and treat 
them as constants in our calculation^. 

The contributions to the polarization tensor of transitions from the nucleon to 
excited states is: 



N' 



Ffj,i/ — Id x 6 



(n\rj N (x)\N) 



(N\r, N (0)\n) ( 



+ 



9(x ) [ (n\ VN (x)\N} § ^ 

+9(x ) (n\ m (x)\N) s ^ =0 [ (N\f} N (o)\n) s ^ Q - (N\f} N (o)\n) s>F ^ =0 ] - 
-e(-x ) [(flMo)!^ - (n\fj N (o)\N) § ^ =0 \ (N\ VN (x)\n)^ =0 - 

-e(-x ) (n\fj N (0)\N) s ^ v=0 [ (N\r)N(xM SiF ^ - (N\r] N (xme, FtiV=0 } (60) 
After collecting terms for each independent tensor basis of the polarization tensor, 



we obtain ( factoring out an overall constant tensor , 
9%\{v) 



2E N M N cos{*f+°f 



9nI(p) 
9nI(p) 
9n\(p) 
9nI(p) 

9m(p) 



><Pn 



r 

On- 



2E N M N sin(^f + ^) 



-2i£#cos(^-^ 
n \ 2 2 



2i [ E^ sin' 
2* ^ 



0n - 
2 ' 

(Vn , On-, 

cos — — H 

M N v 2 2' 



pB-ffN N 
- L N Sm {^ 



2i 



0, 



sin 

M N v 2 2 ' 



E N 



(61) 
(62) 
(63) 
(64) 
(65) 
(66) 



9 A more detailed parameterization including the nucleon pion continuum, which is of importance 
in terms of both chiral symmetry [18] and quantitative error analysis [19], can be used to check 
the validity of this approximation [20]. 
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(c) transitions among excited states 

We shall apply the duality ansatz to identify the contributions to the polarization 
tensor of transitions among excited states; the leading contributions from the 
quark-gluon representation, starting from the continuum threshold Sq , are given 

by 

Re 9% (P 2 ) ( continuum ) = 1 Pr. f I™ g§l(s) ( t^k-gluon) ^ 

The procedure is similar to the case of the nucleon propagator. The quark-gluon 
calculation will be discussed in the next section. 

With all these ingredients, the hadronic representation of the polarization tensor 
Htf(p) is given by the sum of Eqs. (0) ~ (|57D, (|61|) ~ fl66|) , and fl67|), where we 
have two nucleon variables, FJ? , ( or equivalently, F N , o>n ) and four excited state 
unknowns, E^, <Px, E^, tp^ which will appear in the phenomenological side of the QCD 
sum rules. 



IV. STUDY OF THE NUCLEON CORRELATION FUNCTION ( NCF ) FROM 
THE QUARK-GLUON DEGREES OF FREEDOM 

1. Basics 

To represent the NCF in terms of the quark-gluon parameters, we use the method 



of operator product expansion ( OPE ) to calculate the NCF in QCD | 21fl . Because 
the NCF is the vacuum expectation value of composite quark operators, the expansion 
series of the NCF is a short distance expansion in the coordinate space, or, a Aq CD /Q 2 
expansion in momentum spaceQ. We truncate the OPE series at dimension six, and 
the Wilson coefficients are calculated, using diagrammatic rules, in the first order of 
quark mass m q and strong coupling constant g s . 

A new feature associated with our problem is the presence of CP violating operators 
in the OPE series of the NCF. For each CP conserving operator, there exists a parity 
partner which coming from a 75 insertion or a dual transformation, e.g., (qq) g q g 

vs. (qj^q) dq q g and {G 2 )o q e G vs. {GG) 0q Q G . This is also true in the case of induced 
condensates, e.g., {qG^ u q) e 9g vs. {qG^^g g . These CP violating condensates 
introduce new unknowns in the microscopic representation of the NCF. Therefore, 
it is important that we know how to evaluate these new condensates such that our 
calculation has predictive power. 



10 This is in contrast to the use of OPE in the deep inelastic scattering ( DIS ) of lepton off nucleon 
target. Where because of the kinematic ( Bjorken limit ) and the state of interest ( nucleon with a 
four-momentum ) , the expansion series of the correlator in the DIS calculation is a light-cone 
expansion. 
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We can classify the CP violating condensates into two classes: the first class con- 
sists of Ua(1) chiral covariant operators, e.g., i(qj5q) 6 g , whose values depend on 
the representation of the QCD Lagrangian; the second class consists of ?7a(1) chiral 
invariant operators, e.g., (GG) e Bg and {q^^qq^^l^q) e 6 . The classification is useful 
because there is a general relation which connects the parity partners in the first class 
condensates. In particular, one can prove the generalized chiral circle theorem: 

Theorem 2 

] 2 r 1 2 

(qTf(Gfj, u )q) e + i(qTf(G flv ) , y5q) g is invariant under U ~a{1) chiral rotations 



Compared with the simplest case in Theorem 1, here V is an arbitrary Dirac matrix 
and f{G^ v ) stands for any gauge invariant function of the gluon field tensor. 



Use of the polar representation for quark condensates, see Eqs.(|36| ) , (|37|) allows the 
combination of (qq) e Ba + (qi5q) gq g G j5 to be written as 

(qq)e 9 ,e G + (^Q)e q ,e G ^ = -Rfi 4 *™ (68) 

The phase convention for 6 G is chosen such that (qq) e da is negative for 9 q = 9 G = 0. 

The generalized chiral circle theorem can be applied to induced condensates of parity 
partners, e.g.: 

QsqG^q, 9,<rt*GiD><l\ ( 69 ) 

g s qo m G a v q -(p,++u), g s q^a m GZq -(ji++u). (70) 

The susceptibility constants associated with these induced condensates are defined as 
follows: 

g s (qG^q) 0q0G = K q F^{qq) eqeG -ik q F^(q lb q) dqfjG (71) 
9s{qisG^q) eqdG = -iZ, q F^{qq) eqfiG + £ q F^{qj 5 q) Mq (72) 
g s qo m G a v q -(jjt<->v)= V q F^(qq)e q ,e G - iVq F »v{qis(l)e q ,e G ( 73 ) 

Using the generalized chiral circle theorem, we have 

Kq = £q, iq = Kg, Vq = Vq- (74) 

If we assume that all the susceptibility constants are proportional to the quark charge 
e q times a flavor independent constant, then we can rewrite the susceptibility constants 

K q = e q K } £ q = e q £, rj q = e q r]. (75) 

As to the second class condensates, we are only able to use the anomalous Ward iden- 
tity to obtain the value of (GG) e 6a , which is the lowest dimensional chirally invariant 
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condensate in our calculationF]. There are many unknown dimension six chirally in- 
variant condensates, e.g., (qj^qq^j^q) e dc , which cannot be related to their parity 
partners using the chiral circle theorem, and we shall follow the general practice to 
factorize these four-quark condensates into products of dimension three quark conden- 
sates. 

Having explained the subtlety of the calculations of NCF in QCD, we can present the 
result in terms of various coefficients of the invariant tensorsQ: 

2. Quark-Gluon Representation of the Nucleon Propagator U^\p) 

The fi{p 2 ) p sum rule receives the contributions shown in Fig.l . 

Fig.l (a) represents the contribution of the operator /; Fig.l (b) represents the contri- 
bution of the operator m q R q cos 9 q ; Fig.l (c) represents the contribution of the operator 

R 2 q ~ (qi^q qYq)e q ,e G - 

The Feynman diagrams contribute to /f (p 2 )e^ P ^ p2 ^ 75 sum rule are given in Fig.2. 

Fig.2 (a) represents the contribution of the operator m q e ldql5 \ Fig.2 (b) represents the 
contribution of the operator R q e %dal5 ; 

3. Quark-Gluon Representation of the Polarization Tensor U^ u (p) 

Out of the six basis tensors in the polarization tensor, only three of them are useful. 
This is because the (1 , 75) and the (p^j 11 — p v ^) ■ 75 sum rules contain infrared 
singularities which indicate operator mixings with induced condensates whose values 
are unknown^ [pTOfl . Consequently, we only list three sum rules below: 

The e fiual3 p a '-ff3 ■ 75 sum rule contains the nucleon magnetic moments and the OPE 
diagrams are shown in Fig. 3. 

Fig. 3 (a) represents the contribution of the operator F^; Fig. 3 (b) represents the 
contribution of the operator m q R q cos OgF^; Fig.3 (c) represents the contribution of 
the operator R 2 q . 

The p(p ,J ' r y u — p Ur y^) sum rule, being a chirally covariant tensor, contains both the 
anomalous magnetic moment F% and the electric dipole moment F^ , or equivalently, 
Fn and a at. The diagrams in the OPE series are shown in Fig.4. 



n Taking the vacuum expectation value of the anomalous Ward identity, we obtain (GG)^ da oc 



m q R q s'mO q . See ^22 



12 We do not list the complete diagrams which contribute to the OPE series of the NCFs. Instead, 
only terms which survive after further simplifications ( see discussions in the next section ) are 
retained in this section. 

13 For example, the tensor operators m q F^ u , rriqF^ and qcr^q, which appear in the cr^ (1,75) 
sum rule, and the operators F^ U GG and qj^qqjuj^q in the (p^~/ u — p u j^) • 75 sum rule, are mixed 
under renormalization group. 
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Fig. 4 (a) represents the contribution of the operator m q e %6qlh F^, Fig. 4 (b) represents 
the contribution of the operator R q e G75 F Al!/ , Fig.4 (c) represents the contribution of 
the operator (qa^q)^^. 

V. ANALYSIS OF QCD SUM RULES ( QSR ) FOR THE NUCLEON 
CORRELATION FUNCTION: PART ONE 

The QCD sum rules ( QSR ) for the NCFs in the presence of an external electromagnetic 
field will be analyzed in this section. The main emphasis will be on 

• the consistency of the QCD sum rules for NEDMs with symmetry constraints on strong 
CP violations @, and 

• an order-of-magnitude estimate of hadronic variables ( e.g., nucleon masses M N , 
residues A^, EM moments /i^, ) extracted from the QCD sum rules. 

This is a primary analysis in the sense that we shall treat the approximate sum rules ( due to 
truncations of the OPE series and the hadronic parameterizations of the spectral functions 
) as identities, and solve for the hadronic observables as unknown variables in the sum rule 
equations. A more complete study will be presented in the next section. 

1. QCD sum rules of the NCFs in the momentum space 

Here we summarize the results of our calculations in the last two sections by listing 
all the invariant coefficient functions of various independent basis tensors^ for the 
proton correlation function in momentum space. The equalities relating hadronic pa- 
rameterizations and the OPE calculations for various invariant coefficient functions are 
referred to as QCD sum rules. It should be kept in mind that these two representations 
of NCFs are derived from different expansions and truncations ( a hadronic complete 
set in the former and 1/Q 2 in the latter ) of the same correlation functions. Therefore, 
these equalities are at best approximate identities whose validity are empirical. 

The QCD sum rules for the proton propagator n£ ' representation are given by 

nf = f?( P 2 )P + /1VK W)7S (76) 

(a) fi{p 2 ) ( p sum rule ) 

K> . p 4 \n(—p 2 ) Am u a u cos 9 u \n(— p 2 ) 2a 2 
1- - + continuum = y — \ g ' + " " K ' - t^t^TT ( 77 ) 



p 2 -M* 4(2tt) 4 (2tt) 4 3(2tt)V 

(b) / 2 p (p 2 )e^( p2 ^ 5 ( J, 7s sum rule ) 



14 We do not list the (7^(1,^) sum rules because of the infrared problem | 
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___ e ^ + continuum = 4(27r)3 e + ^ e ^ (78) 

The QCD sum rules for the proton polarization tensor n^,(p) in the polar form rep- 
resentation are given by 

n^) = E^V)-T* (79) 



(c) #!<y ) ( ^ uaP PaJ(3 ■ 7s sum rule ) 



2AW + F|) _ 2A?f^cos(4 - t) 



4(p 2 - M 2 ) 2 



+ 



4(p 2 - M 2 ) 



= e d 



p 2 ln(-p 2 ) + m u a u cos fl M + a 



2(27if 



3p 4 



+ continuum 



+ (oc e u ) 



(80) 



(d) #4 (p ) ( (p M 7^ — p^7 M ) • 75 sum rule ) 

2A 2 Ff 2A 2 [^sin(4-|)-^sin(f -I)] 

H - — — — — r h continuum 



4(p 2 -M 2 ) 2 4(p 2 -M 2 ) 
(oc m q R q sin 9 q ) (81) 



(e) #f(P ) + Wl{p )ls ( p(p"7" - PV) sum rule ) 



-2)A 2 F P 



4M p (p 2 — M 2 ) 5 



3 i(9 p +a p )75 



+ 



4M p (p 2 - M 2 



)75 



67. 



m d ln(-p 2 ) 
(2^) 4 



+ 



Xad hi(-p 2 ) a d 
6d 3(2vr) 4 + 6u (2vr) V 



+ continuum 



(82) 



Here we define a g = (27r) 2 -R, « 0.55(Ge\/) 3 . 

The neutron sum rules can be obtained from the proton ones by doing an isospin 
rotation, namely, replacing the hadronic observables, e.g., M p by M n , X p by A n ; and 
the QCD parameters, e.g., by m u ,and Rd by R u etc. 

2. The use of the Borel transform 

In principle, the QCD sum rules, as summarized above, can be used to extract infor- 
mation on hadronic variables in terms of QCD parameters. However, the validity of 
the matching between these two representations for the NCFs is severely limited by the 
convergence property of the OPE series and the uncertainty of the higher state contri- 
butions. Fortunately, there are several prescriptions to improve the sum rule relations. 
Specifically, we wish to generate a set of improved sum rules such that the OPE se- 
ries have better convergence and the contributions from a hadronic parametrization 
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is dominated by the ground state observables. For this purpose, we use the Borel 



transformation, which is defined as [16 

B MB [f(Q 2 )]^ 



lim 



Q 2 -> oo 
n — > oo 
M B = ^ fixed 



-HQ 

n 



2\n+l 



d 



n+1 



dQ 2 



(83) 



We need to apply this transformation^ on both sides of the sum rule relations: 
For instance, on the OPE side of the QCD sum rules, we have 

B M J(Q 2 ) m \n(Q 2 )} = (-l) m+1 m!(M| 
fel 1 




(*-l)! 



Ml, 



(84) 
(85) 



Thus, after Borel transformation, the higher dimensional operators in the OPE series 
receive further suppression with factorial factors. 

On the phenomenological side of the QCD sum rules, we have 



B 



M B 



Q 2 + M 2 



. M 



{k-l)\ 



(86) 



The power suppressions of the excited state contributions are Borel transformed into 
exponential ones. 

Having established the usefulness of the Borel transformation, we need to specify how 
the matching of sum rule relations can be realized. As we have mentioned before, 
the identities we have derived cannot be exact for all values of M|. A choice of 
"matching region" has to be made; such a choice is a compromise between different 
convergence properties of the two representations of the NCFs. On the OPE side, 
we prefer a large value of Mb to suppress power corrections ( see Eqs.(^), (|85"D ); 
on the phenomenological side, we prefer a small value of Mb to enhance ground state 
observables ( see Eq . (|86"D ). In view of this, the matching of the sum rule relations only 
works for a finite range of the Borel mass M#, and the region of matching the QCD 
calculations and hadronic parametrizations is generally referred to as a Borel window. 

It is an empirical fact that within the Borel window, physical quantities have a mild 
dependence on the squared Borel mass M B . Hence, in the primary analysis, we can 
choose a value for the squared Borel mass M B close to the ground state mass Mjf to 
obtained a rough estimation for the observables of interest. 



15 Notice that, after the Borel transformation, the virtual 4-momentum variable Q 2 
pearing in the momentum space sum rules is replaced by the squared Borel mass M B . 



-p 2 ap- 
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A technical remark regarding the continuum contribution: 

The continuum contributions to both nucleon propagator and the polarization tensor 
are given by the duality ansatz Eq. (|46|) , fl67|). After applying the Borel transform 
Eq.(||D ( with k = 1 ) on both sides of Eqs.fff), we get 



B 



M B 



R- e fi (P )( Continuum ) 



Pr 



Im f i (s)( quark-gluon 



'b ds 



(87) 



The leading OPE terms of the NCF are of the form of (p 2 ) n ln(— p 2 — ie), which has an 
imaginary part equal to — ir(p 2 ) n . Thus the right hand side of the Eq.(|87|) is a truncated 
Laplace transformation of a polynomial in p 2 . Since the continuum contributions share 
a similar form as the leading OPE terms on the quark-gluon side of the NCF, we can 
combine these two terms together, and define the following functions: 



E n (s,w) = 

E (s,w) = 
E^s, w) = 

E 2 (s,w) = 



I - / dt e" t n 

w 



e 
e 



(88) 

(89) 
(90) 

(91) 



All functions E n act as a high-energy cutoff for the leading OPE series, with E n (s, w) w 
if s > w. Since the continuum threshold represent an average parameter for the 
excited state spectrum, we shall choose w = Sq = 1.75GeV 2 and write E n {M^) = 
E n (s = M 2 B) w = 1.75GeV 2 ). 



3. Borel transform— improved sum rules 

The Borel transform-improved sum rules for the proton propagator IIo(p) is given by 
(a) /f (Ml) 



M%E 2 {M 2 B ) + Am u a u cos 9 U M 2 B E (M 2 B ) + 



Aa 2 



X 2 p e 



(92) 



(b) f 2 (M 2 B )e^ v{ v ^ 



2M B E 2 {M 2 B )m d e w ^ + 2M B E 1 (M 2 )a u e 



i(>Gdl5 



X 2 p M p e ^ e ^ 75 



(93) 



where = 2(2n) A X 2 y . For the proton polarization tensor 11^ (p), we obtain the 
following sum rules: 

(c) 9 P 3 (M 2 ) 
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e u M^Ei(M 2 B ) + e d m u a u cos 9 d + 



4[-(e d + f) + (oce tt)] 



A 



n + n 



Ml 



+ E£ cos 



0„ 



(94) 



[ 4e u a d + (oc e d ) \ < 
F, 



5 - [ Ae u m d + (oc e d ) 



M 2 B e 



Ml 



P e i(d p+ap )-y 5 + E B e i(^l^L) l5 



(95) 



4. Manipulation of the QCD Sum Rules 

While we have identified the hadronic unknowns and written down the sum rule equa- 
tions to relate hadronic observables in terms of quark-gluon parameters, we need to 
take some further simplifications to make the solutions of these complicated relations 
manageable^. 

The three steps we adopt to simplify the sum rule relations are discussed in order: 

(a) The use of isospin symmetry 

It is important that the symmetry constraints on the strong CP violations require 
that all quarks have finite masses. In particular, the existence of any massless 
quark leads to no strong CP violation. In such a limit, the flavor symmetry is 
broken and we have to make the flavor dependence in our calculation explicit. 
While it is possible to realize these limits in the QCD sum rule approach ||22|| . 



we opt for isospin symmetry to simplify our calculations. Consequently, for the 
quark-gluon parameters, we take m u = m d = m q and R u = R d = R q . It can be 
shown that in this case, 8 U can be chosen to be equal to 9 d and, as a consequence 
of Crewther's condition, 6 Gu = 6 Gd = 6 G /2. See Eqs.([3Bp ([59]). 

On the other hand, since we do not include QED corrections, which break the 
isospin symmetry, in the calculations of the hadronic masses and residues, we can 
take M p = M n and \ p = X n . Thus, the use of isospin symmetry greatly reduces 
the unknowns and the input parameters in the sum rule relations. 
The other advantage of this simplification, as pointed out by B.L. Ioffe and A.V. 



Smilga |i~0 |) is that we can eliminate the induced condensates which generate 
the unknown susceptibility constants from the proton and neutron sum rules by 
taking certain combinations with different quark charge dependences. 



"However, we emphasize that none of these simplifications are absolutely necessary for the phys- 
ical content of the further discussions. The choice of these simplifications should be regarded as 
pedagogic and to help us isolate issues of strong CP violation from other complications in the 
strong interactions. 
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For example, in the p{p^^ v — p u, y fl )(l, 75) sum rule we can multiply the pro- 
ton sum rule Eq.(p2|) by e u and subtract the corresponding neutron sum rule 
multiplied by e^. This eliminates the contributions of the induced condensates 
(Wwdoifia ^ x{qq)e g ft G - Fo r the eP^pcflp ■ 75 sum rule Eq.(f3), we repeat the 
similar procedure, but with times the proton sum rule minus e u times the 
neutron one. As a bonus, we also eliminate the infrared singularity coming from 
the operator rriqF^. 

(b) The elimination of excited state parameters in the polarization tensor 

We shall not concern ourselves with the details of the excited states in this calcu- 
lation. With the assumption that the ground state to excited state transition can 
be approximated by momentum independent constants Eq.^, one can observe 
that such a contribution has a different Borel mass dependence, as compared to 
the ground state observables, see Eqs. Q94]) (P5p. We can apply the differential 
operator 







1 - (96) 



on both sides of the sum rule equations Eqs. ([54]), (j9o| ) to eliminate the contribu- 
tions of ground to excited state transitions in the sum rule relations. 

The problem of operator mixing in the {p^ r ) v — p 1/ 7 /i ) • 75 sum rule 
After all these simplifications, we have seven unknowns A at ,9 n ,F p ,a p ,F n ,a n in 
the sum rules relations. While the nucleon propagator gives three indepen- 
dent sum rules ( p, 1, 75 ) for the nucleon variables A at, Mjy, 6^, we have only three 
useful identities Eqs. (|9^) , (p5|) in the polarization tensor 11^, for four EM moments 
F p , a p , F n , a n \ The leading OPE contribution to the tensor (p^j 11 — p Ur y^) ■ 75 is 
proportional to the operator F flu (GG) q &g , whose Wilson coefficient has a lnm q 
infrared singularity^. One can argue that such an infrared singularity should be 
included in the definition of QCD condensates, which requires a separation scale 
instead of a small quark mass. The difficulty here is to show our calculation is 
independent of the choice of a separation scale, as required by the renormalization 
group analysis. The situation is further complicated by the mixing of operators 
under renormalization. Without a detailed analysis of the operator mixing prob- 
lem, we can not obtain useful information from this tensor. Consequently, we 
need to introduce a phenomeno logical parameter (3 for the ratio of F% to F£. By 
comparing the Feynman diagrams contributing to the tensor (p fl '-f u — p U/ -f tl ) ■ 75 
for both proton and neutron, we find the only difference is the charge dependence 
of u and d quarks, the quark masses and chiral radii do not enter this sum rule 
because of the chiral property of this sum rule. The charge dependence seems 
to indicate that (3 = e u /ed- Henceforth, we shall assume such a relation in our 
subsequent discussion. 
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We need to calculate a quark propagator in the presence of two or three constant external fields 
( G^y and/or F^ u ), which is infrared divergent if the quark mass is zero. 
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5. Symmetry Constraints from the QCD Sum Rule Relations 

Despite the complicated appearance of the sum rule relations, it is possible to show 
that, due to the use of the polar form in both hadronic and quark-gluon representa- 
tions of the NCF, our sum rule calculations satisfy the symmetry constraints on the 
strong CP violation B without an explicit solution of the CP violating observable, 
e.g.,NEDM. 

To see this, we focus on the (1, 75) sum rules, Eq.fl9~3|) in the nucleon propagator and 
the p(p flr j u — p v j fl )(l, 75) sum rules, Eq.(p5[) in the polarization tensor. Both sets of 
sum rules are associated with chirally covariant tensors and receive contributions from 
chirally covariant condensates. In a shorthand notation, we can rewrite the sum rules 
in the following forms: 

A m q e w ^ 5 +B R q e i(>G ^ 5 = Ce idN ^ 5 (97) 
A' m q e^ 75 + B' R q e ie °^ = cV (ejv+ajvhs (98) 

We emphasize that the structure of the sum rules, as summarized above, is independent 
of the approximation scheme which only affects the values of the numerical coefficients 
A,B,A f and B' . Specifically, in terms of the Wilson coefficients appearing in the ( 
Borel transformed ) QCD sum rules, these numerical coefficients are 



A = 2M B E 2 (M 2 B ) (99) 

B = 2M B E l (M 2 B ) (100) 

A' = 4(4 - el)M B E (M 2 B ) (101) 

B' = -4(e 2 , - el)Ml (102) 



In the two sets of sum rules Eqs.(|97|) (|98"D, we have four hadronic unknowns C,C',9n, 
and «jv to be determined. The relation with the nucleon variables is that C is a 
function of and Mjy, C is a function of F p and F n , with an assumption^ that 

tan" 1 -! = tan" 1 -! = a N (103) 

p n 

we rearrange the seven nucleon variables into four unknowns. 

The former sum rule, Eq.(|97|) contains the nucleon chiral phase 8n, the latter Eq.(|9~8"D 
has both 9^ and a^, in the hadronic representations of the NCF. It is the relative ( 
chiral ) phase difference that defines the physically measurable ( and chirally invariant 
) dimensionless ratio d n / [i a w which determines the violation of the CP symmetry. 

On the other hand, the OPE series for both sum rules are organized in such a way 
that the explicit chiral symmetry breaking parameter m q and the spontaneous chiral 



8 This assumption is chosen for its simplicity and is not equivalent to the conjectured value for 

/3 = e u /e d . 
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symmetry breaking parameter R q are on equal footings^. Such a representation of 
the NCF is particularly useful for examining the symmetry constraints on the strong 
CP violation. Namely, in either limit m q — > or R q — > 0, the nucleon chiral phase 
9 N becomes 9a q or 9 q , respectively, and the relative phase = tan _1 (c? ri //x^) has 
to vanish. Thus, there is no strong CP violation if chiral symmetry is exact in QCD 
and CP violating observables must be proportional to the product of m q and R q . 
Furthermore, the periodic structure in the chiral phases naturally generates a sva.9 
factor^, which combined with the two chiral symmetry breaking parameters m q and 
R q , implies that CP violating observables are proportional to (GG)^ 

We can solve the equations fl9~T|), (|98| ) for the hadronic unknowns C,C',9N,aN as 
follows: 

C 2 N = A 2 m\ + B 2 R\ + 2ABm q R q cos 9 (104) 

C% = A' 2 m 2 q + B' 2 R 2 + 2A'B'm q R q cos 9 (105) 

Am q sin 9 q + BR q sin 9 G 

t&n9 N = — — -p — (106) 

Am q cos 9 q + BR q cos 9q 

t n = (AB' - BA')m q R q sm9 

an a N AA , m 2 + BB'R\ + ( AB' + BA')m q R q cos 9 1 J 

The nucleon variables are given by the following expressions: 

G 

m _^ = r 108) 

N M 6 E 2 {M 2 ) + 4m q R q M 2 E {M 2 ) cos # g + 4B 2 /3 1 J 

A^v = M B E 2 (M B ) + 4m q R q M B cos9 q E0(M B ) + 4/^/3 (109) 
M 2 

= ^4(e„f/(Ml,^) - e d K(M B ,9 q )+e 2 ) (110) 



M 2 

F 2 = -r-^iedUiMl, 9 q ) - e u K(M 2 B , 9 q ) + e d e u ) (111) 



M 



Fi = -^_e u V(M Bl 9 q ) (H2) 
M 2 

F 3 n = —^e d V(M 2 B ,9 q ) (113) 



where 



9 The form of the chirally covariant sum rules is a concrete realization of the strong CP torus of 
QCD, as described in refer. ||. 

20 This is the simplest periodic function, which is regular and approaches zero as 6 vanishes. 

21 Taking the expectation value of the anomalous Ward identity for the flavor singlet axial current, 
one can show that (GG)g is propotional to the product of m q , R q and sin 9. 
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K(M 2 B ,9 q 



U(M 2 B ,9 q 



V(M 2 B ,9 q 



I -Ml 



d 



dM B 



B dM 2 



B . 



B . 



I 2 2 



AR 2 q /3 + Am q R q M 2 E (M 2 ) cos9 q 



M 6 E 2 (M B ) + 4m q R q M 2 E {M 2 ) cos9 q + AR^ q /3_ 
AA'm 2 + (AB' + A'B)m q R q cos 9 q + BB'R 2 



N 



{AB>-A'B)(GG) dqfiG 



2C% 



(114) 
(115) 
(116) 



Given the small value of the current quark mass, m q « 5MeV, and the relative minus 
sign in the numerator of Eq. (|107|) , it is possible that the CP violating observables 
could be dynamically suppressed without a tiny 9 parameter. We shall examine such 
an interesting scenario in the next section. 

6. Order of Magnitude Estimations of the Hadronic Observables from the 
QCD Sum Rules 

Before we embark on a complete analysis of the sum rule relations, it is useful to 
estimate the relative sizes of the hadronic observables. The qualitative features of the 
estimations also provide a test whether the underlying assumptions in the sum rule 
calculation work for the case in which we are interested. It is an empirical fact that 
the Borel window lies around the lGeV region and we can choose a value for the 
Borel mass Mb to be the nucleon mass ~ lGeV, if the Borel mass dependence of the 
hadronic observables is weak. Other inputs needed for the estimations of the hadronic 
observables include the current quark mass m q w 5MeV, and the quark chiral radius 
R q « -(qq)e q =e G =o = (2A0MeV) 3 . The latter gives a q = 0.55(GeV) 3 . 

Substituting these numbers into Eqs.(^) ~ ( |100| ), in the small 9 limit, the functions 
Eqs.( TO ~ ( PD reduce to 



M N {9) 
\ 2 {9) 
F p 2 {9) 
F?(9) 
Fi(9) 

mo) 



0.9GeV 
2.5(GeV) 6 
3.1 
-2.1 



10 -2 sin 9 



-0.5 x l(T 2 sin# 



(117) 
(118) 
(119) 
(120) 
(121) 
(122) 



Comparing the ratio /F^ 1 with the current experiment upper bound on the nEDM 
d n / — 10~ u , we obtain an upper bound on the strong CP violating parameter 
9 < 10~ 9 . 

We emphasize that the upper bound on the 9 parameter is obtained from a calculation 
without assuming a perturbative expansion on the 9 parameter. Consequently, we 
can derive an upper bound on the ratio d n / [i^ with respect to the 9 parameter. Such 
information is an intrinsic property of the CP conserving (9 = 0) QCD and provides 
a dynamical suppression mechanism to the CP violating observables for solving the 
strong CP problem. Unfortunately, the number we have derived max 0< g <27r % < 
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10 2 , is not small enough to achieve this end. Therefore, we conclude that the solution 
to the strong CP problem has to lie beyond QCD. 

VI. ANALYSIS OF QCD SUM RULES ( QSR ) FOR THE NUCLEON 
CORRELATION FUNCTION: PART TWO 

This section aims at a more careful study of the QSRs for the NCFs in the presence of an 
external EM field. We shall extract various nucleon variables, e.g., the nucleon masses Mat, 
nucleon residues and their EM moments , F^ 1 , from the QSRs in a more rigorous 
manner. The emphasis here is to use the QCD sum rule method to obtain quantitive 
results for the hadronic variables, including an analysis and/or estimate of the errors and 
uncertainties in our calculations. 

We first discuss how to choose a region for the matching between two representations of 
the NCFs; then we extract the hadronic observables from a least square fit. The possible 
errors and uncertainty are discussed in the third part and we summarize the results obtained 
from the sum rule analysis in the last part. 

1. Choice of the Borel Window 

As we have emphasized before, the hadronic and quark-gluon representations of the 
NCFs are based on two different expansion schemes and there is no a priori reason that, 
under our truncation and approximation, these two expansions should provide the same 
information of the NCFs. While it is possible to enlarge the overlap and improve the 
convergence with the use of the Borel transform, a matching region between both sides 
of the sum rules has to be specified. In view of the different convergence properties, 
we use the following criterion to define the Borel window. 

First of all, on the quark gluon representation of the NCFs, we require that the highest 
dimension condensates contribute to the whole series less than 20%. This leads to a 
lower bound on the Borel mass squared M\ > 0.8 GeV 2 . 

Secondly, for the hadronic parametrization, we restrict the contribution from the con- 
tinuum to be no greater than 20% of the leading OPE contributions to the NCF. This 
leads to an upper bound on the Borel mass squared M| < 1.2 GeV 2 . 

Thirdly, the extractions^ of the hadronic variables in the region specified above intro- 
duce additional uncertainity. Presumably, the third uncertainity is not independent 
from the ones associated with the convergence requirement. However, such a correla- 
tion is not transparent in our analysis, and we simply take the intersection region of the 
three criterions to define the Borel window, which lies between 0.8GeV 2 to 1.2GeV 2 . 

2. Extraction of the Hadronic Observables from the QCD Sum Rules 



There is no standard procedure to extract hadronic variables in the sum rules approach. In our 
case, we average the hadronic variables over the Borel window and use the x 2 to represent the 
uncertainity with our extractions. 
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Given the simplifications discussed in the previous section, we obtain several relations 
which give the hadronic observables as functions of the squared Borel mass. As we have 
emphasized, if the matching between the hadronic and quark-gluon representations 
of the NCF is realized in the given Borel window, all the hadronic observables should 
have only a mild dependence on the squared Borel mass. Therefore, we plotQ all six 
chirally invariant hadronic observables, nucleon mass M/v, nucleon residue X%, proton 
and neutron magnetic moments 1 + F$, F% and electric dipole moments F%, F£, as 
functions of the squared Borel mass Mj, in Fig. 5 and Fig. 6. The final values of these 
hadronic observables are determined by averaging the functions over the Borel window 
and the uncertainty is given by the x 2 °f this average. Since it is necessary to have 
such x 2 uncertainty to be smaller than the uncertainty associated with the choice of 
Borel window, we shall consider the result obtained from our sum rule analysis to lie 
within 20% uncertainty. 

3. Error Analysis 

We list a few comments on the errors and uncertainties related to this calculation. 
On the OPE calculations of the NCF: 

(1) The convergence of the OPE expansion: 

The OPE series for the NCF in momentum space is a Aq CD /Q 2 expansion, where we 
assume that all vacuum condensates are proportional to Aq C d to a certain power, with 
a numerical coefficient of order 0(1). While there is no systematic way to estimate 
the large order behavior of the coefficients and test this assumption, the expansion 
series could be an asymptotic series and has no convergence radius at all. We have to 
content ourselves with the naive estimate from the Borel window analysis and estimate 
the contributions from the higher dimensional condensates to be below 20%. 

(2) Approximations in the calculations of the Wilson coefficients: 

When we calculate the OPE series of a correlation function, the Wilson coefficients are 
expanded in the powers of "small" parameters, at s , m 2 q /Q 2 and we make truncations 
of the series to obtain an approximate results for the Wilson coefficients. Again, the 
higher order behavior of the perturbative expansion is poorly known and we make no 
attempt to study the higher order contributions in the present work. 

(3) Hadronic parameterizations: 

This is the most challenging part in the sum rule analysis, because we have no ad hoc 
criterion to see how well our ansatz, e.g., the momentum independence of the excited 
state parameters E% , Eg, works. There are some analyses which show that this part 
could bring in a large uncertainty to a sum rule calculation [fL9| . 

(4) The extraction of physical quantities from the matching: 

Since physical quantities should be independent of the Borel mass, the deviation of the 
physical observables from a constant value in a given Borel window should not be larger 



In presenting these figures, we have chosen the value of the 9 parameter to be 10 
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than the uncertainty associated with the choice of a Borel window. Nevertheless, the 
scheme dependent errors ( how do we extract physical quantities from the matching 
) are related to the x 2 of the matching, and represent the quality of a sum rule 
calculation. 

(5) The dependence on the input parameters: 

The QCD sum rule calculations generally rely on input parameters which are derived 
from other sources, e.g., the quark mass, quark condensates etc.. In many cases, these 
numbers are not well determined and their variations from a "standard" value have 
to be taken into account when we estimate the errors of the physical quantities from 
a QCD sum rule calculation. Recently, an analysis was performed by D.B. Leinweber 
[[2~3|j and we shall not dwell on this aspect further in this work. 

4. Summary of the Results 

In summary, we obtain the following results for the hadronic observables from the sum 
rule analysis: 

M N = 0.9GeV (123) 
X 2 N = 2.5GeV 6 (124) 
Fl = 3.1 F™ = -2.1 (125) 

if = 1(T 2 F 3 n = -0.5 x 1(T 2 (126) 



VII. SUMMARY AND CONCLUSION 

The electric dipole moments of nucleons ( NEDM, ) are calculated using the method 
of QCD sum rules. Through the use of a polar representation of both nucleon EM moments ( 
anomalous magnetic moment and electric dipole moment F^ ) and the {7a (1) covariant 
quark condensates, we are able to demonstrate the reparameterization invariance of CP 
violating observables explicitly in a QCD Lagrangian with a CP violating 9 parameter. The 
symmetry constraints on strong CP violation in QCD, together with a dual relationship 
between the quark mass m q and the chiral radius R q |§ are realized in the QCD sum rule 
relations which connect the hadronic observables to the QCD parameters. The extraction 
of the NEDM in terms of the 9 parameter and QCD parameters from the QCD sum rule 
relations can be achieved without assuming a perturbative expansion of the 9 parameter. In 
addition, the periodic dependence of the CP violating observables on the 9 parameter comes 
out naturally in our approach. 

Our final result establishes a functional dependence of the NEDM on the 9 parameter, 
which, combined with the experimental upper bound on the nEDM, gives an upper bound on 
the unknown 9 parameter of less than 10~ 9 . While this result is compatible with previous 
calculations of the nEDM using different techniques, it also indicates that a dynamical 
suppression of the CP violating observables, as implied by the symmetry constraints on the 
strong CP violation in QCD, is not sufficient to resolve the strong CP problem. Therefore, 
we conclude that a solution to the strong CP problem does not exist within QCD and a 
natural explanation of the tiny 9 parameter necessarily requires physics beyond the standard 
model. 
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FIGURES 

FIG. 1. Feynman Diagrams which contribute to the p sum rule 
FIG 2. Feynman Diagrams which contribute to the 1 and 75 sum rules 
FIG 3. Feynman Diagrams which contribute to the t^^Pctlfilh sum rme 
FIG. 4. Feynman Diagrams which contribute to the p(p fJ, j u — p u j^)(l, 75) sum rule 
FIG. 5. The nucleon mass Mn as a function of the squared Borel mass M\ 
FIG. 6. The nucleon residue as a function of the squared Borel mass M\ 
FIG. 7. The nucleon magnetic moment F® + F% as a function of the squared Borel mass M\ 
FIG. 8. The nucleon electric dipole moment F^ as a function of the squared Borel mass 
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TABLES 

TABLE I. Nucleoli to higher state EM transition matrix element 
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